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Hadronic scattering mediated through the strong interaction has been an area of great interest for 
both theory and experiment. Recently, lattice QCD calculations of scattering processes from first 
Q_i principles have seen remarkable progress, including multiple calculations of two-pion scattering 

lengths that agree with experiment to within a few percent. However, there exists a certain class of 
scattering processes, such as pion-nucleon scattering, that contain annihilation diagrams, which 
are often prohibitively expensive to simulate on the lattice. In this talk, I will present a method 
to extract certain parameters from this class of scattering processes by employing an isospin 
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determinant. 
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1. Introduction 

Within the past decade, there has been substantial progress in understanding scattering phe- 
nomena from first-priciple lattice QCD calculations through analyzing two-hadron correlation func- 
tions. Several meson-meson calculations, most notably the 1=2 %% scattering channel, agree with 
experimental calculations to within a few percent [1, 2, 3, 4]. In addition to the meson-meson 
channel, the multi-baryon channel has received a great deal of attention [5] and continues to show 
considerable progress. 

One set of processes that has not received as much attention are baryon-meson processes. 
These processes are needed not only to understand long-range interactions of nucleons through 
pion exchange, but kaon-nucleon interactions play a pivotal role in understanding potential kaon 
condensation in neutron stars [6]. While calculations have been preformed for several meson- 
nucleon scattering processes [7], there exist a class of processes that have eluded these calculations; 
namely, calculations with annihilation diagrams. While immense progress has been made in recent 
years due to advancements in algorithms and the use of GPUs, calculation of annihilation diagrams 
are often prohibitively expensive requiring all-to-all propagators. 

In this work, we present an alternative method of extracting these parameters through the use 
of an isospin chemical potential. This method, proposed in Ref. [8], takes advantage of the dy- 
namically generated pion condensate when the isospin chemical potential reaches a certain critical 
value. Then, through the aid of heavy baryon effective field theory in its region of applicability, the 
relevant scattering parameters can be extracted by performing baryon spectroscopy in the presence 
of this pion condensate. 

2. Annihilation Contributions in Meson-Baryon Processes 

As mentioned previously, currently lattice methods of extracting scattering information rely on 
two (or more) hadron correlation functions. For meson-baryon processes, this consists of putting a 
meson and baryon operator in the source and sink, performing the relevant contractions, and then 
summing over all possible final states. For some meson-baryon processes, such as 7T + £ + , there are 
only "non-annihilation" processes, as shown in Fig. 1. As a result, only a point-to-all propagation 
(effectively one column of the matrix) is necessary to extract the relevant information, and these 
methods were preformed in Ref. [7] . 

However, there exist another class of scattering processes, such as tcN or K~N, that contain an- 
nihilation processes, which can have a significant effect on the scattering process (the A-resonance 
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Figure 1: Examples of the "non-annihilation" contributions to the 7T + E + process. 
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Figure 2: Examples of the "annihilation" contributions to the n + n process. 

in K~N scattering is an example of this). These effects, as illustrated in Fig. 2, require an all-to- 
all propagation (whole matrix is required), which is often very expensive to calculate using brute 
force methods. In light of this issue, we explore a different calculation, that only requires baryon 
spectroscopy in the chiral regime. 



3. Taking advantage of condensates 

It has been well established for many years that condensates can be used to shift or generate 
masses. The Higgs mechanism in particle physics is an example of this process as the Higgs 
vacuum expectation value gives masses to the vector bosons. So could this process help extract 
scattering processes on the lattice? For simplicity, we explore a system where a scalar without 
flavor indices interacts with a baryon. The term in the Lagrangian for this process is given by 

jSf = aB<p 2 B+--. (3.1) 

In most low-energy processes that we are interested in, the value of a is a function of the non- 
perterbative physics that can only be extracted through lattice calculations or experiment. As men- 
tioned previously, this process may contain annihilation diagrams that may render the lattice cal- 
culations prohibitively expensive. However, if we were to dynamically form a scalar condensate, 
namely (j) (x) = 0o + (x), the term in the Lagrangian becomes 

Jz? = a^BB + ■ ■ ■ . (3.2) 

Since a<j>Q is a constant, this term can be thought of as a shift of the baryon mass, as shown in 
Fig. 3. To be more precise, if the baryon mass without the condensate is M^-q then the baryon 

mass in the presence of the condensate would be M^ )= q + a<pQ -\ . If we knew and could control 

the value of 0n> we could, in principle, extract the parameter a simply by taking differences of the 
baryon masses. 

So the two questions that remain are, first, do we have a mechanism to dynamically generate 
a condensate, and, two, if we generate this condensate, do we have a method to systematically 
account for its effects? Indeed, there are methods that allow both of these questions to be answered 
favorably. The isospin chemical potential can generate a pion condensate and effective field theory 
gives us a controlled tool to see how the mass will shift as a function of this condensate. 
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Figure 3: Diagrammatic illustration of how the condensate can shift the mass. By taking the difference of 
the masses calculated in both pictures and varying 0o, one can extract the scattering parameter of interest. 

4. Isospin chemical potential 

One well known problem that occurs when studying systems at finite baryon density is the 
fermion sign problem. When the baryon chemical potential is non-zero, the Euclidean action no 
longer has a positive definite determinant, which prohibitively hinders the usual Monte Carlo meth- 
ods for lattice simulations. As a result, there has been much effort studying finite density systems 
without a fermion sign problem. One such example is the isospin chemical potential. 

The action for the isospin chemical potential for the two flavor theory in Minkowski space is 
given by 



By inspection (after the Euclidean rotation), one can see that the usual condition of 75-Hermiticity is 
not satisfied. However, a different condition Ti/s-hermaticity 1 is satisfied, which ensures a positive 
definite determinant: 



This feature has enabled multiple lattice calculations at finite isospin density [9, 10, 11]. It is 
also worth pointing out the symmetries of this action. When m = jJ. = 0, the action processes the 
usual SU (2)l <g> SU (2) # flavor symmetry which spontaneously breaks to SU (2) y . When m/0 and 
;U 7^ 0, the action processes a U(1)l i S)U(1)r flavor symmetry which breaks, both explicitly and 
spontaneously to £/(l)y. 

In addition to not having a prohibitive fermion sign problem, the isospin chemical potential 
leads to a pion condensation phase [9, 10, 11, 12, 13]. This condensation phase, which breaks 
the U(l)v symmetry spontaneously, occurs when the chemical potential reaches a certain critical 
value. Basic arguments imply that this transition would occur when pL =m n , since there are no 
lighter 1=1 particles that could be generated for smaller chemical potential values. However, for 
a theory where the quark mass is light enough for chiral perturbation theory to be valid, this can be 
shown explicitly. 

1 T 2 75-Hermiticity is equally valid for this Lagrangian. 




(4.1) 





(4.2) 
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5. Connections to effective field theory 

The way you map the isospin chemical potential term to the chiral effective field theory is to 
treat the chemical potential term as a gauge field, namely 

T 3 

V^My^o, (5.1) 
while obeying the power counting p 2 ~ m\ ~ \i 2 . The resulting chiral Lagrangian is given by 

f 



if 



tr(D M J/D' i J/ t ) + 2Atr(M t £/ + C/ f M) 



(5.2) 



where D^U = d^U + i[Y^,U] and U is an SU (2) matrix of the pseudo-Goldstone bosons and can 
be written as 

U = e iaAz = cosa + /ir-Tsina. (5.3) 
Minimizing the potential yields the result [12, 13] 

2 

m 

cosa = ^. (5.4) 

There are several concepts that can be taken away from this equation. First, when jj. < m K , the 
system dynamically behaves the same way as the /i = scenario and when jj. > m n , a pion con- 
densate is formed and the chiral condensate gets altered. Second, both the chiral condensate and 
pion condensate can be given in terms of the parameter a. If we add a small isospin breaking term 
to ensure that the spontaneous breaking of the pion condensate is in the T2 direction (a process that 
is conveyed in more detail in Ref. [8]), the values of the condensates are given by 

(xj/y) = / 2 Acosa, (5.5) 
/(y7r 2 y 5 i//) = / 2 Asina. (5.6) 

Throughout the rest of this note, the results will be given in terms of cos a and sin a. However, 
thanks to these relations, we can relate cos a to the value of the chiral condensate and sin a to the 
value of the pion condensate. 

To reiterate the goal of this work, we want to find a relation between the mass shift of a hadron 
in the presence of an isospin chemical potential to the value of the condensate, chemical potential, 
and scattering parameters. Thus, if we know the value of the condensate and chemical potential, 
along with the mass shift calculated on the lattice, we should be able to extract several elusive 
scattering parameters. It should also be reiterated that the measurements of these mass shifts do not 
require all-to-all propagation and can be calculated using the usual spectroscopy methods. In order 
to find this analytic relation between baryon mass shifts and the condensate/chemical potential, 
heavy baryon chiral perturbation theory is employed [14]. It should be emphasized at this point 
that for this analytic relation to be valid, the effective field theory must not break down. In other 
words, these relations will only be useful for pion masses and chemical potential values in the 
chiral regime. 
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We will now present the results and implications for the two flavor nucleon calculation 2 . The 
nucleon mass dependence on the chemical potential and condensates is given by (in terms of the 
coefficients from the heavy baryon lagrangian in Ref. [14, 15]) 

M N = M^ - jU/cosa^- +4ci (mucosa + Xe sin a) + ( c 2 - §fj + c3 ) jU/ 2 sin 2 a, (5.7) 
2 V oM J 

where the low energy constant c\ is the "sigma term" and c 2 and C3 determine the scattering length 
for pion-nucleon scattering. The e in this equation represents the small isospin breaking term 
needed to point the pion condensate in a particular direction in flavor space. As mentioned previ- 
ously, this is explained in Ref. [8]. The key message to take away from this equation is that the 
chemical potential (and consequently, the condensates) acts as an additional "knob" for extracting 
low-energy constants. Starting with the sigma term, c\, which is still difficult to determine pre- 
cisely on the lattice [17], by varying jU, the value of cos a will vary, and c\ can be extracted at a 
fixed pion mass. This is a particularly interesting result as you can extract information about chiral 
physics if the pion mass were at the physical point. The scattering parameters c 2 and C3 can also 
be extracted by plotting the mass shift of the nucleon vs the chemical potential times the pion con- 
densate. It should be pointed out, that this method alone cannot disentangle c 2 and C3. However, 
the implementation of twisted boundary conditions could separate these terms. Additionally, it is 
worth mentioning the approximate size of the mass shift and whether or not it can be reasonably 
extracted. If we use the estimate of ci ~ .9 GeV- { [15] it is not difficult to see the nucleon mass will 
need to be determined from .5% to 5% depending on the pion mass (smaller pion masses require 
better resolution in the spectroscopy calculation). 

6. Conclusions and Future Applications 

The primary motivation for exploring the effect of isospin chemical potentials on hardronic 
mass shifts was to extract information about scattering processes that require calculation of annihi- 
lation diagrams on the lattice. This all-to-all process is very expensive and this particular method 
allows one to move the issues from the propagator calculations to generating multiple configura- 
tions at different chemical potentials (in a similar way lattice calculations are generated at multiple 
mass values). However, this approach can be more appropriately though of as an additional "knob" 
on the theory that one can extract chiral information from. In Ref. [8], the baryon-meson channel 
was the primary focus, but a similar calculation could be preformed for the meson-meson chan- 
nel where calculating meson masses in the presence of an isospin chemical potential will lead to 
meson-meson scattering. 

This begs the question of whether or not other potentially interesting information can be ex- 
tracted from lattices with an isospin chemical potential. Since pion condensation is one of the 
possible phases that could exist within neutron stars, it would be worthwhile to understand the 
low-energy phenomena in this state. Some potentially interesting calculations include heavy quark 
potentials, two meson phase shifts, etc. Also, it would be enlightening to explore the thermody- 
namics of these systems. If one can vary the temperature and isospin chemical potential, it should 

2 For the complete analysis, along with the calculations including the strange quark, see Ref. [8], using the two flavor 
hyperon analysis in Ref. [16] and performing the relevant matching. 
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be possible to find a region of small fugacity when the virial expansion is valid. The coefficients in 
these expansions contain information about three and four meson interactions and could be directly 
extracted from the lattice using these methods. 

References 

[1] T. Yamazaki et al. [CP-PACS Collaboration], "I = 2 Pi Pi Scattering Phase Shift with Two Flavors of 
0(A) Improved Dynamical Quarks," Phys. Rev. D 70 (2004) 074513 [arXiv:hep-lat/0402025]. 

[2] S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage [NPLQCD Collaboration], "I = 2 Pi Pi 
Scattering from Fully-Dynamical Mixed- Action Lattice QCD," Phys. Rev. D 73 (2006) 054503 
[arXiv:hep-lat/0506013]. 

[3] S. R. Beane, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, A. Torok and A. Walker-Loud, "Precise 
Determination of the 1=2 Pipi Scattering Length from Mixed-Action Lattice QCD," Phys. Rev. D 77 
(2008) 014505 [arXiv:0706.3026 [hep-lat]]. 

[4] X. Feng, K. Jansen and D. B. Renner, "The Pi+ Pi+ Scattering Length from Maximally Twisted Mass 
Lattice QCD," Phys. Lett. B 684 (2010) 268 [arXiv:0909.3255 [hep-lat]]. 

[5] S. R. Beane et al. [NPLQCD Collaboration], "High Statistics Analysis Using Anisotropic Clover 
Lattices: (Iii) Baryon-Baryon Interactions," Phys. Rev. D 81 (2010) 054505 [arXiv:09 12.4243 
[hep-lat]]. 

[6] D. B. Kaplan and A. E. Nelson, "Strange Goings on in Dense Nucleonic Matter," Phys. Lett. B 175 
(1986) 57. 

[7] A. Torok et al, "Meson-Baryon Scattering Lengths from Mixed- Action Lattice QCD," Phys. Rev. D 
81 (2010) 074506 [arXiv:0907.1913 [hep-lat]]. 

[8] P. F. Bedaque, M. I. Buchoff and B. C. Tiburzi, "Meson-Baryon Scattering Parameters from Lattice 
QCD with an Isospin Chemical Potential," Phys. Rev. D 80 (2009) 1 14501 [arXiv:09 10.4595 
[hep-lat]]. 

[9] J. B. Kogut and D. K. Sinclair, "Lattice QCD at Finite Isospin Density at Zero and Finite 
Temperature," Phys. Rev. D 66 (2002) 034505 [arXiv:hep-lat/0202028]. 

[10] J. B. Kogut and D. K. Sinclair, "The Finite Temperature Transition for 2-Flavor Lattice QCD at Finite 
Isospin Density," Phys. Rev. D 70 (2004) 094501 [arXiv:hep-lat/0407027]. 

[11] P. de Forcrand, M. A. Stephanov and U. Wenger, "On the Phase Diagram of QCD at Finite Isospin 
Density," PoS LAT2007 (2007) 237 [arXiv:071 1.0023 [hep-lat]]. 

[12] D. T. Son and M. A. Stephanov, "QCD at Finite Isospin Density," Phys. Rev. Lett. 86 (2001) 592 
[arXiv:hep-ph/0005225]. 

[13] D. T. Son and M. A. Stephanov, "QCD at Finite Isospin Density: from Pion to Quark Antiquark 

Condensation," Phys. Atom. Nucl. 64 (2001) 834 [Yad. Fiz. 64 (2001) 899] [arXiv:hep-ph/00 11365]. 

[14] V. Bernard, N. Kaiser, J. Kambor and U. G. Meissner, "Chiral Structure of the Nucleon," Nucl. Phys. 
B 388 (1992) 315. 

[15] M. Frink and U. G. Meissner, JHEP 0407, 028 (2004) [arXiv:hep-lat/0404018]. 

[16] B. C. Tiburzi and A. Walker-Loud, Phys. Lett. B 669, 246 (2008) [arXiv:0808.0482 [nucl-th]]. 

[17] A. Walker-Loud et al, Phys. Rev. D 79, 054502 (2009) [arXiv: 0806 .4549 [hep-lat]]. 



7 



